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Abstract: The effects of parton energy loss in cold nuclear matter on heavy-quarkonium 
suppression in p-A collisions are studied. It is shown from first principles that at large 
quarkonium energy E and small production angle in the nucleus rest frame, the medium- 
induced energy loss scales as E. Using this result, a phenomenological model depending 
on a single free parameter is able to reproduce J/V' and T suppression data in a broad Xp- 
range and at various center-of-mass energies. These results strongly support energy loss 
as the dominant effect in heavy-quarkonium suppression in p-A collisions. Predictions for 
i /tp and T suppression in p-Pb collisions at the LHC are made. It is argued that parton 
energy loss scaling as E should generally apply to hadron production in p-A collisions, 
such as light hadron or open charm production. 
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1 Introduction 

The spectacular quenching of hadrons produced at large in Pb-Pb collisions at the 
LHC [1, 2], as well as the jet imbalance reported in those collisions [3, 4], find a natural 
explanation in terms of parton energy loss in a quark-gluon plasma (QGP). For light hadron 
production at mid-rapidity and sufficiently large , the parton energy loss is dominantly 
radiative, in average of the form AE ~ ghot [5, 6], with L the distance travelled by the 
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parton through the hot medium and ghot the rate per unit length of transverse momentum 
broadening in the medium. The strength of jet-quenching can be explained if the transport 
coefficient g in a QGP is larger, by one to two orders of magnitude, than its estimate in 
cold nuclear matter, t^coid ~ 0.045 GeV^/fm [7]. This is why jet-quenching is considered as 
a prominent QGP signal. However, despite the wealth of data accumulated so far at RHIC 
and LHC, the in-depth understanding of energy loss processes in a QGP remains far from 
complete (see [8] for a discussion). 

Drastic nuclear suppression effects are not only seen in A-A but also in p-A collisions, 
at least for some processes and in some kinematical conditions. For instance, quarko- 
nium [9] but also light hadron [10, 11] production at large longitudinal momentum fraction 
Xp (or large rapidity) is strongly suppressed in p-A as compared to p-p collisions. Under- 
standing nuclear suppression in cold nuclear matter, a well-controlled medium as opposed 
to an expanding QGP, should be a prerequisite in order to interpret quantitatively nu- 
clear suppression in heavy-ion collisions. However, it is striking that there is no consensus 
yet on the origin of i/ip suppression at large rapidity/x^ in p~A collisions, from SPS to 
RHIC [9, 12, 13], despite many theoretical attempts (see [14] for a review). 

Recently, new scaling properties have been identified for the induced gluon radiation 
spectrum dl/dcj, and associated energy loss A£', of hard processes where a color charge 
undergoes small angle scattering through a static medium (cold matter or QGP) [15]. In 
the present work we address the phenomenological consequences of these results on i/ij) 
and T nuclear suppression in p-A and vr-A collisions, parametrized by the attenuation 
factor (in the following we use the generic notations "-0" and "p-A") 

< = ^ ^ ^) • (1-1) 

We will show that the available large-x^ quarkonium suppression data in p-A collisions 
can be explained by parton energy loss in cold matter. Although q in cold matter is small, 
a strong nuclear attenuation arises due to the specific parametric behaviour (x E at 
sufficiently large E, where E is the quarkonium energy in the target nucleus rest frame. As 
discussed in Ref. [15] and reviewed in the present paper (Section 2), this behaviour holds 
when the hard partonic subprocess can be viewed, in the nucleus rest frame, as the small 
angle scattering of a color charge. In the following we focus on quarkonium hadroproduction 
(see Fig. la), where the heavy quark mass provides the hard scale allowing for a perturbative 
QCD description, and for which p-A suppression data are quite abundant. Our discussion 
should however apply more generally to hadron hadroproduction, for instance to light 
hadron production in p-A collisions (provided the light hadron p± plays the role of the 
hard scale, i.e., p± ^ 1 GeV), see Fig. lb. Light hadron nuclear suppression due to parton 
energy loss will be addressed in a future work. 

As is well-known, the quarkonium hadroproduction mechanism in elementary p-p col- 
lisions is still under debate^ (see for instance Ref. [18] for a review). In order to study 

^This is especially true for low p± production, p± < M, considered in the present study. The recent 
findings of Refs. [16, 17] apply specificaUy to the large p± 2> M domain. 
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Figure 1. Generic processes of (a) heavy-quarkonium liadroproduction (b) light liadron liadropro- 
duction (c) deep inelastic scattering and (d) DrcU-Yan production, at large E in the target nucleus 
rest frame. The ellipse represents the hard subprocess occurring within the time thard- Cases (a) 
and (b) are similar to small angle scattering of an asymptotic charge. 



quarkonium nuclear suppression in the most model-independent way, in the present study 
we will only assume that the heavy-quark QQ pair of mass M is produced, within the per- 
turbative proper time scale Tqq ~ 1/M, in a compact color octet state, and remains color 
octet for a time Toctet ^ ''"qq- quarkonium production models where color neutralization 
is a soft, non-perturbative process, Toctet coincides with the quarkonium hadronization time 
r^, and this assumption holds at any Xp. In the Color Singlet Model (CSM), the gluon 
emission required for color neutralization of the QQ pair is constrained at large enough 
Xp to become softish and thus to occur late, leading to Tqq <C Toctet ^ t^- Thus, the 
assumption of a color octet QQ pair living longer than the perturbative time scale ^ 1/M 
holds quite independently of the quarkonium production model. ^ 

Working in the nucleus rest frame and considering the limit E ^ M ^ p±, quarko- 
nium hadroproduction looks like small angle scattering of an "asymptotic" color charge, 
i.e., prepared in the "far past" and propagating in the "far future" as compared to the 
perturbative time scale thard = Tqq " (E/M) ~ E/M"^. This is illustrated in Fig. la for the 
generic gg — )■ QQ partonic subprocess, viewed in the nucleus rest frame as the splitting 
g —7- QQ of the incoming gluon, followed by a rescattering in the target.^ 

In the present study we will assume for simplicity that the octet QQ pair arises dom- 
inantly from the splitting of an incoming gluon. This should be a valid assumption for all 
p-A data considered in this paper, except at very large values of Xp {xp ^ 0.8), where 
quark-induced processes come into play. (We will further comment on this point in Section 

^Strictly speaking, the assumption is independent of the quarkonium production model only at large 
Xp. At small Xp, the assumption becomes invalid in the CSM. We will briefly come back to this point in 
our final discussion. 

^If the gluon in Fig. la is relatively large, gx ~ C {M), then it should be interpreted as part of the 
hard subprocess gg — >■ QQ + g (or gq — >• QQ + q) . 



-3- 



4.1.) 

Within those assumptions, the associated gluon radiation with large formation time, 
tf ^ ihard) is similar to (non-abelian) Bethe-Heitler radiation off a fast color octet under- 
going an effective transverse momentum kick q±. The typical q± is expected to be larger in 
p-A than in p-p collisions due to transverse momentum broadening being proportional to 
the target size L, Aq'j_ = ~ qL.^ As a result, the medium-induced radiation spectrum is 
similar to the Bethe-Heitler spectrum (2.8), up to the replacement of the total momentum 
transfer by the broadening through the nucleus A (see (2.17) and Section 2), 

with M± = (M^ +P^)2 the transverse mass of the QQ pair and Ap = max(A^^j^,£^p). 
This leads to an average loss AE oc E. When A^ < £^ . ^ M? we have 



AE^ / duju—c^^N.as^^-^ E. (1.3) 

Jo dcj M_L 

The spectrum (1.2) is at the basis of the phenomenological study presented here, 
whose main results can already be found in Ref. [19] . The scaling AE (x E in quarkonium 
hadroproduction was first postulated in [20] (also revisited in [21]) yet this assumption was 
not motivated and the parametric dependence on L and M arbitrary (and different from 
(1.3)). In Ref. [22], an energy-independent bound on AE was derived, but in a specific 
setup where the nuclear broadening of the final tagged particle was neglected (we will 
briefly comment on this point in Section 2). 

We stress that the spectrum (1.2) is coherent. This can easily be seen in a calculation 
using physical polarizations for the radiated gluon (see Section 2). With this choice the 
medium-induced radiation spectrum indeed arises from the interference between the initial 
and final state emission amplitudes. The energy loss (1.3) is thus neither a purely initial 
nor final state effect, and is distinct from gluon radiation resummed in leading- twist parton 
distribution and fragmentation functions. Being process-dependent {e.g., it is expected in 
3/ip production in p-A collisions but not in Drell-Yan production, see below) and sup- 
pressed by a power of the hard scale M±, it is naturally interpreted as a higher- twist effect. 
Nevertheless it plays a crucial role in a broad or rapidity interval, as we shall see. 

To make the physics under consideration clear, let us mention that the spectrum (1.2) 
is not expected in quarkonium (real or virtual) photoproduction, nor in inclusive deep 
inelastic scattering (DIS) off nuclei, where the incoming energetic particle participating 
to the hard subprocess is colorless (see Fig. Ic for the DIS case).^ In those cases radia- 
tion with tf ^ ^hard can only arise from final state radiation. The latter (DGLAP-like) 
radiation is independent of the medium properties and cancels in the medium-induced 



*In our study the broadening Aq^ equals the amount of soft rescattering In other words it is defined 
with respect to an "ideal" target where soft rescatterings are absent, see Section 2.2. 

^This is to be distinguished from resolved photoproduction, which at the partonic level is similar to 
hadroproduction (Fig. la), and where we thus expect AE oc E. 
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spectrum. Quite remarkably, at large z the quarkonium production data in deep inelastic 
muon scattering [23] exhibit no nuclear suppression (but instead a slight enhancement), in 
sharp contrast to hadroproduction. 

Drell-Yan (DY) production off nuclei is similar to quarkonium photoproduction, since 
by definition the energetic particle produced perturbatively (the Drell-Yan photon of mass 
Q) is colorless, see Fig. Id. Radiation with if » thard ~ E/Q"^ must arise from initial state 
and does not contribute to the medium-induced spectrum. Thus, neither the spectrum 
(1.2), nor the energy loss A£J oc E' of the type (1.3), is expected in DY production. This 
is qualitatively consistent with the much milder nuclear suppression of DY production [24] 
when compared to i/ip hadroproduction in the same kinematical range. 

In Section 2 we discuss the medium-induced gluon radiation spectrum, derived in [15], 
which is at the basis of the model detailed in Section 3. Phenomenological applications are 
presented in Section 4 and we conclude the paper by a discussion (Section 5). 



2 Revisiting energy loss scaling properties 

In this section we justify the expression (1.2) for the gluon radiation spectrum associated 
to quarkonium hadroproduction, in a less heuristic way than in Ref. [15]. 

2.1 Asymptotic charge 

We first consider the case of an on-shell ("asymptotic") parton of energy E undergoing 
an elastic scattering and exchanging a gluon with transverse momentum (.±_ with a nuclear 
target, see Fig. 2a. As is well-known, scattering can induce radiation, provided the quantum 
state of the charge is perturbed. In QED, this happens when the scattering angle Os — 
is non-zero. In QCD, radiation can occur even in the limit Og — )• 0, due to the incoming 
parton color rotation in the elastic scattering. 



E 




(iii) 

(a) (b) 

Figure 2. (a) Elastic scattering amplitude M.^i of an asymptotic light quark via single gluon 
exchange, (b) Induced gluon radiation amplitude AArad- 

These features are illustrated by the expression of the gluon radiation amplitude in- 
duced by the elastic scattering, given by the diagrams of Fig. 2b. We denote by u and 
fc_L the radiated gluon energy and transverse momentum. We focus on soft {u <C E) and 
small angle (fcj^ ^ <^) gluon radiation. For an on-shell light quark the radiation amplitude 
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reads [25] 
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(2.1) 



(2.2) 



with e the physical polarization vector of the radiated gluon, which will be implicit in the 
following. The first two terms of (2.1) correspond to initial state radiation (diagrams (i) 
and (ii) of Fig. 2b), and the last (diagram (iii) of Fig. 2b) to final state radiation. 

In the abelian case, the second diagram is absent, as well as color factors, and the 
radiation spectrum reads 
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(2.3) 



The collinear divergence at — )• is screened by the electron mass m. Keeping the latter, 
the radiation spectrum can be expressed as^ 
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(2.4) 



where 6m = m/E. The radiation spectrum vanishes when i± = 0, but also in the formal 
limit 0s — )• (at fixed Om), as expected. 

In QCD we single out the purely non-abelian contribution to the radiation spectrum 
off a light quark by focussing on the 0^ — )■ limit of Eq. (2.1), 
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(2.5) 



Squaring the r.h.s. of (2.5), we recover the well-known Gunion-Bertsch spectrum [25], 

di e 



duj d^/c I 



2.6) 



Integrating over d'^k±, or equivalently over cPO, and inserting the quark mass M dependence 
as in the abelian case above, we obtain^ 



dl_ 

duj 



Nr.a, 



QCD 



d02 



(2.7) 



^Although the mass dependence of the QED spectrum (2.4) is not exact, it provides the correct para- 
metric expressions in the two limits £± ^ m and £± m [15, 26]. This remark also applies to the QCD 
expression (2.7). 

''The color factor in (2.7) is obtained by summing \Mrad\^ over initial and final color indices, and 
normalizing by the color factor associated to |7V/fe;i^- The resulting factor Nc is independent of the type 
(quark or gluon) of the energetic color charge. 
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where 9m = M/E and the lower cutoff arises from the constraint fc_L > A (with A = Aq^^), 
put by hand for the consistency of our perturbative QCD treatment (we also assumed 
i±_ > A). We obtain the soft radiation spectrum off an on-shell quark of mass M [15, 26], 



The average energy loss of a heavy quark, M ^ £±, is dominated by cj ~ {i±/M) E <^ E 
and is proportional to E, 

rE J T / _ A 

AE= duju— ~ N,a,^-—E. (2.9) 

Jo d'^ A/>f± M 

2.2 Color charge resolved in a hard process 

The case of a color charge resolved in a hard process can be simply illustrated by inserting 
a hard exchange q± » £± in the light quark scattering process, see Fig. 3a. The transfer i± 
now plays the role of nuclear momentum broadening, £^ = Ag^. Since some radiation is 
released even when i± = due to the presence of the hard exchange, the relevant quantity 
is the medium-induced radiation spectrum. 



d/ 

dcj 



^ uj^{q±; i±) - uJ^iqr, = 0) . (2.10) 
ind du 




(a) (b) 

Figure 3. (a) Model for hard scattering q± supplemented by transverse momentum broadening 
£j_ <C q±. (b) Associated gluon radiation amplitude in the large formation time and soft gluon 
limit. 

We concentrate on soft radiation as compared to the hard process, oj <^ E, k± <^ 
and on the limit of large gluon formation time. It can indeed be checked that the domain 
if ^ max(L,thard) contributes most to the medium-induced loss (see Section 2.3). The 
radiation amplitude is thus dominated by the diagrams of Fig. 3b. Diagrams of the type 
shown in Fig. 4 can be neglected. The diagram of Fig. 4a, relevant when tf < L, contributes 
to the energy loss as AE oc L^, which is negligible compared to the contribution AE oc E 
we focus on throughout our study. The diagram of Fig. 4b is suppressed due to i± <^ q_L. 

Note that in the limit tf ^ L, the induced radiation cannot probe the nuclear size L, 
but is however sensitive to the total amount of soft rescattering, = qL, transferred to 
the fast color charge, independently of the actual number of soft rescatterings. Modeling 
soft rescatterings by a single effective (semi-hard) scattering with £^ = qL should thus 
allow addressing the main features of medium-induced radiation. 
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Figure 4. Diagrams for tlic radiation amplitude which arc suppressed at large formation time and 
in the approximation k±,£± ^ q±. 





(a) 



(b) 



Figure 5. Interference term (a) and associated virtual correction (b) contributing to the medium- 
induced spectrum (2.15). 



The radiation amplitude is similar to the amphtude (2.1) in absence of hard scattering, 
however with modified color factors and scattering angle Og, 
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(2.12) 



In the abelian case, the second term of (2.11) is absent. Neglecting l± compared to 
q± in the definition of Og, the radiation amplitude becomes independent of and the 
induced spectrum (2.10) vanishes. As argued by Brodsky and Hoyer [22], this means that 
only radiation with small formation time tf < L can contribute, resulting in some {E- 
independent) bound on medium-induced energy loss. This statement holds in an abelian 
model (such as that considered in Ref. [22]), but fails in QCD, as we shortly recall now.® 

In QCD the radiation spectrum is obtained by squaring (2.11) and putting aside purely 
initial/final state radiation, which cancels in the medium-induced spectrum. The latter thus 
arises from the interference between the second and third diagrams of Fig. 3b, see Fig. 5a, 
and reads 
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(2.13) 



In the following we maintain the approximation Os\^f^ — ^s|pp = l^/E, which allows to extract the 
purely non-abelian medium-induced spectrum (1.2). This approximation corresponds to an experimental 
setup where the p-p and p-A distributions are measured in the same px-bin. If the distributions 
are averaged over px, we expect Os\^p^ > 6s|pp, leading to some abelian-hke radiation arising from U 2> L. 
However, the associated loss is suppressed by a power of the hard scale as compared to the purely non- 
abelian loss. 
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The azimuthal integral yields [15] 
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ind Jdg " Jmax(2:2g-[,A2) 



dl 










do; 


ind ^ 





where x = uj/E and we added the constraint k± > A in the A;_|_-integral. Approximating 
{xq± - i±)'^ ~ x'^q'j_ + i]_ and niax{x'^q'j_, A^) ~ x'^q'j_ + A^ we obtain 

It is easy to check that for the scattering of a fast octet charge (the case of interest), the 
color factor in (2.15) reads = Nc- Note that the virtual correction shown in Fig. 5b 
contributes a factor 2 in (2.15).^ 

We can also verify that with a parton mass M ^ we get the same expression as 
(2.15), up to the change in the hard scale qj_ — )■ M± = (M^ + q'j_)^ [15]. The spectrum off 
an energetic color octet charge thus reads 



d/ 

do; 



ind ^ L \ ^"_L 



?1 



where the 0-function reminds us that only ^_|_ > A can induce the emission of perturbative 
gluons. Quite remarkably, the medium-induced spectrum (2.16) is parametrically similar 
to the radiation spectrum (2.8) of an asymptotic parton of "mass" M±. 

Finally, we stress that in the above calculation, the medium-induced spectrum has 
been defined with respect to an "ideal" target for which £± = 0, see (2.10). In practice, 
nuclear suppression is measured in a nucleus A with respect to that in a nucleus B < A, 
where soft rescatterings can also occur (even in the proton case B = 1). The medium- 
induced spectrum relevant to this situation should be defined with respect to the nucleus 
B and reads (as already announced in (1.2) for B = 1) 

(...) 

where A| = max(A2^^,.^^g). When the broadening in the nucleus B is too small, ^5_b ^ 
^QCD' spectrum (2.17) coincides with the spectrum (2.16) defined with respect to the 
'idea 

B becomes independent of A, 



"ideal" target with £± = 0. When > A^ , the induced radiation in A with respect to 



QCD • 

2.3 Application to quarkonium hadroproduction 

As illustrated by the previous section, the behaviour AE oc E for medium-induced parton 
energy loss is not forbidden by first principles. We expect such a behaviour in all p-A 
processes where some color charge is scattered at small angle in the nucleus rest fame, in 
particular in quarkonium hadroproduction at large E. 



^For a quick way to derive color factors and the relative contributions of virtual corrections, see for 
instance Ref. [27]. 
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However, in order to apply the spectrum (1.2) to quarkonium hadroproduction, the 
underlying partonic subprocess should effectively look like the scattering of an energetic, 
pointlike color charge, at least within the formation time if ^ L of the radiated gluon. 
This is the case if 

max(L,thard) < if < ioctct < ^^"^ r_i_{tf) <^ l/k± , (2.18) 

where thard is the hard process time scale, toctet the lifetime of the color octet QQ pair, 

and the quarkonium hadronization time (see the Introduction). The second condition 

states that the QQ pair is effectively pointlike when the transverse wavelength l/k± of the 

radiated gluon is larger than the transverse size r± of the quark pair at a time ~ tf. In 

principle, the two conditions (2.18) can be checked a posteriori, using the typical formation 

time contributing to the observable of interest (in our case nuclear attenuation). 

As an illustration, we show that the typical tf contributing to the average loss (1.3) 

formally satisfies (2.18). For simplicity we assume here L < thard ~ E/M'j_ and Ad^ ~ M, 

and denote l\ ~ It is trivial to check that in (1.3) arises from radiated energies 

(jj r-j E{i±/M) <C E and transverse momenta k"^ ~ (see (2.14)). The typical tf thus 

satisfies „ „ ^ 

E Lo E E 

the inequalities arising from the nuclear broadening i± being soft compared to M, and hard 
(for large nuclei) compared to the non-perturbative scale t^^ ~ 0.6 GeV (see Section 3.5.1). 
In this respect, let us recall that momentum broadening is related to the saturation scale Qs 
in the nucleus, = Q^ [28], indeed considered as a semi-hard scale. The second condition 
in (2.18) reads 

k± r^(tf) ~ k^ tf ~ £^ • ■ ~ a. « 1 , (2.20) 

where the relative transverse velocity of the heavy quarks is estimated by v± ~ Pb±/E, 
with pb ~ OsM the Bohr momentum of the quarkonium state. Thus, the conditions 
(2.18) are fulfilled in the perturbative domain ^ 1 and provided nuclear broadening 
is a semi-hard scale. This defines the theoretical limit where the energy loss (1.3) can be 
applied to quarkonium hadroproduction. 

In practice, £± ~ VqL » t^^ might not be satisfied. Indeed, using q = 0.08GeV^/fm 
(see Section 4) and L = 7fm we have y/qL ~ 0.7 GeV ~ ^ somewhat questioning the 
validity of (1.3) for quarkonium hadroproduction. However, it should be stressed that 
the observable of interest in this paper - nuclear attenuation - does not directly depend 
on the average energy loss, but rather on the energy loss probability distribution V{ijJ,E), 
see Section 3. As is well-known and generic to jet-quenching phenomenology [7], nuclear 
attenuation is dominated by the low energy tail of V{uj,E), i.e., by uj much smaller than 
the typical uj contributing to the average loss. This leads to smaller values of tf to be used 
in (2.19), leading to the required condition tf <^ t^. 



^°The QQ pair is produced perturbatively at the time thard with a relative momentum Sp± < M between 
the quark and antiquark. However, only those pairs with 5p± < UsM eventually have a non-negligible 
overlap with the quarkonium wave function. 
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Under the conditions (2.18) we thus expect the induced radiation spectrum to be given 
by Eq. (1.2), as derived in Section 2.2 for a fast color octet charge. It should be clear from 
Section 2.2 that the parametric dependence of the spectrum is uniquely determined, and 
should thus apply to other processes than quarkonium production (like open charm and 
light hadron hadroproduction), as discussed in the Introduction. 

Finally, we note that within our approximation (2.18) the quarkonium bound state 
is formed far beyond the nucleus, 3> L. This approximation may break down at low 
proton beam energy or at small values of Xp. We will comment on this when discussing 
the limits of applicability of the model in Section 3.5.1. 



E,V^)= deV{e,E)^{E + e,^) . (3.1) 



3 Model 

3.1 Shift in energy or "medium-induced splitting" 

The starting point of the model consists in expressing the J/V' differential production 
cross section da/dE in p-A collisions simply as that in p-p collisions, with a shift in 
the quarkonium energy E accounting for the energy loss e incurred by the octet cc pair 
propagating through the nucleus, 

A dE 

The energy loss e is more conveniently defined in the nucleus rest frame, and we thus 
denote E and Ep ~ s/{2mp) the J/V' and projectile proton energies in this frame (with irip 
the proton mass). We have e < Ep — E from energy conservation, and we impose e < E 
for consistency with the soft radiation approximation. Hence emax = min(£^p — E, E) in 
(3.1). The quantity V{£,E) is the energy loss probability distribution or quenching weight 
associated to the radiation spectrum (1.2), to be discussed in Section 3.4. 

The measured differential cross sections are usually expressed as a function of (or 
of the rapidity) rather than E. The variable x^ is defined as the longitudinal momentum 
fraction between the J/ip and projectile proton in the cm. frame of an elementary p~N 
collision (of energy y/s). It can be related to the J/ip transverse mass M± and rapidity y' 
in this frame. In the limit y/s ^ nip, 

_ P'w 2M^ sinhy' 



Using E = M_|_coshy = M±cosh.{y' -\- Ay), where cosh Ay = ^/s/{2mp) with Ay the 
projectile proton rapidity in the cm. frame, we obtain from (3.2) 



E = E{Xp) = Ep 
The relation (3.3) can be inverted to give 




(3.3) 



, , E Er, M'j , , 

(3.4) 
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Introducing the variable x', 

E Ep M'] 2M± 



X 



' - x'{E) = ^ + ^ = cosh y' = Jxl + AMl/s , (3.5) 



E 



we have dx^/dE = x' /E and obtain from (3.1) 

'^\x^,V^)^f-.eVie,E) [(||^] ^ (..(^ + .), V^) • (3.6) 



A dx. 



Note that at large Xp » M^/^/s, we have Xp ~ x' ~ E/Ep, the (Jacobian) factor 
in between brackets in (3.6) approaches unity, and the p-p cross section is evaluated at a 
shifted value of x^, with the shift 6xp{e) = Xp{E + e) — Xp{E) ~ e/Ep. At large Xp the 
energy shift in (3.1) is thus equivalent to a simple translation in x^. This is not true at all 
values of Xp , due to the presence of the Jacobian. 

In the following we will use the expression (3.6), where E = E{Xp), and the relations 
(3.3) and (3.4), valid at all Xp. In a model where the J/V' is produced through a 2 — ?• 1 
partonic subprocess, the expression (3.3), denoted as E = xiEp, simply arises from the 
standard relations between parton momentum fractions, xiX2 = Mj_/s and xi — X2 = 
Xp (note also that xi + X2 = x'). However, the kinematical relation (3.3) is actually 
independent of the partonic subprocess. 

It is interesting to mention that the main equation of our model (3.1) is equivalent to 



where z = E/[E + e) is interpreted as a (medium-induced) splitting variable describing 
the energy loss process {E is the energy of the charge after radiating the energy e), and 
F\oss{z) as a "medium-induced splitting function". The expression (3.7) follows from (3.1) 
by changing variable from e to z (g iving ZYain — max(£'/£/p, 1/2)), and using the fact that 
the quenching weight V{£,E) is a scaling function of the ratio ejE^^ 

e\ ^ /I - 2 



E ■ V(e, E)=r[-)=r l^-^ j=z' Jioss(^) . (3.8) 

In writing (3.1), we implicitly assumed that the energy of the radiating octet cc pair is 
the same as the final J/V' energy E. However, the equivalent expression (3.7) suggests that 
our nuclear suppression model based on a simple shift in E might apply to more general 
situations where the final detected particle's energy arises from the fragmentation of some 
parent parton's energy, with a fragmentation variable z < 1. Indeed, suppose that the 
observable p-p cross section is of the form^^ 

d.DM^K'^V (3,9) 



Je/Ep dE \ z 



^^This can be trivially checked from (3.23) (Section 3.4) and (1.2). 

^^The dependence of (3.9) on ^/s and on the projectile and target parton distribution functions is irrelevant 



to our discussion. 
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and that medium-induced radiation and hadronization factorize/^ 

where z'^^^ = max{E/{zEp),l/2). Changing the order of the z and z' integrals in (3.10) 
and using (3.9), we recover (3.7). 

The quarkonium p-A cross section is thus related to the observable p~p cross section 
according to (3.7), or equivalently (3.1), independently of the form of the fragmentation 
function D{z). This result mostly follows from the scaling property (3.8) of the quenching 
weight, and is thus expected to hold in all processes involving a. fractional medium-induced 
energy loss (A£^ oc E), in particular in open charm and light hadron hadroproduction. 
Instead of describing the energy loss process as a shift in energy (see (3.1) or (3.6)), one 
could alternatively describe it as in (3.7) by a medium-induced splitting function^'^ T\oss{z)- 

Finally, although in the above discussion we assumed the p~p cross section to obey 
factorization (see (3.9)), we believe that (3.1) might hold independently of this assumption. 
For instance, as long as the underlying partonic process is similar to the scattering of 
a color charge (as in Figs, la and lb), we may imagine the p-p cross section at large 
E to be affected by late comover rescattering [30] and thus to violate factorization, and 
nevertheless the p-A cross section to be given by (3.1). The only crucial assumption is that 
the partonic subprocess induces radiation, as dictated by perturbative QCD, independently 
of the precise mechanism fixing the quantum numbers of the final detected particle. 

3.2 Absolute production cross section 

The dynamics of heavy-quarkonium production in hadronic collisions is still uncertain. In 
particular, none of the existing models proposed to describe heavy-quarkonium production 
is able to reproduce simultaneously all the features reported experimentally, at both pj_ < 
M and p± > M. 

In the present approach, a crucial ingredient entering (3.6) is the single differential 
absolute cross section, dapp/dXp, of J/^/J and T production in p-p collisions at a given 
center-of-mass energy. In order to be as model-independent as possible, dapp/dxp used in 
(3.6) is not taken from theory but determined from a fit to the data. We found that it can 
be conveniently parametrised as 

_i;(,,,„L_i, (3.11) 

where x' is defined in (3.5) and the exponent n is extracted from p-p^^ and 7r~-p data 
taken at the same center-of-mass energy - whenever possible - as the p-A and 7r~-A 

^^This should be guaranteed by the separation of time scales, t( <^ t^, see Section 2.3. 

^^This designation is motivated by the fact that Juoss is perturbatively calculable and can be Mellin 
convoluted with the (vacuum) fragmentation function D{z) to give the "medium modified fragmentation 
function" Dmcdiz) — J dz' D{z') Fioaa{z / z') . We however stress that J^ioss(z) is process-dependent {e.g., it 
is present in quarkonium but absent in DY production, see the Introduction). Thus, D^ed[z) in the latter 
equation differs from the medium modified fragmentation functions assumed to be universal and discussed 
elsewhere, see for instance Ref. [29]. 

^^p-A data, where A is a light nucleus such as Be or C, were also used, see Fig. 6. 
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p beam 



TT beam 



Experiment 


NA3 


E789 


HERA-B 


PHENIX 


ALICE 


NA3 


NA3 


NA3 


^/5 (GeV) 


19.4 


38.7 


41.5 


200 


7000 


16.8 


19.4 


22.9 


n 


4.3 


4.5 ±0.05 


5.7 ±0.2 


8.3 ± 1.1 


32.3 ±7.5 


1.4 


1.4 


1.5 



Table 1. Values of n extracted from J/ip production in p-p (or p-Be, p-C) and tt -p collisions; 
see text for details. 



Experiment 


E866 


PHENIX 


LHCb 


(GeV) 


38.7 


200 


7000 


n 


3.4 ±0.2 


6.7 ± 1.0 


14.2 ±2.9 



Table 2. Values of n extracted from T production in p~p (or p-d) collisions; see text for details. 

measurements discussed in this paper. ""^^ Note that the normalization parameter in (3.11), 
or equivalently the total production cross section, is irrelevant for our purpose since only 
nuclear production ratios are considered in this paper, see Eq. (1.1). 

The values of n for J/^/; production from SPS to LHC energy are summarized in 
Table 1.^^ The index for p-p production grows smoothly from n ~ 4-5 at low energies 
(y^i < 40 GeV) up to n ~ 8 at RHIC and n ~ 30 at the LHC. At LHC, the relative 
uncertainty on n is as large as 6n/n ~ 20% {n = 32.3 it 7.5) because of the too small Xp 
domain, \xp \ < 0.02, covered by the data and the fact that around mid-rapidity, x' <^ 1, the 
parametrization (3.11) becomes independent of n. However we checked that the resulting 
uncertainty on i?p^ at LHC (and similarly at RHIC) is marginal because of these very 
reasons. In the T channel, the exponents are slightly smaller than for J/ip production; see 
Table 2. 

The comparison between the fits and J/V' (respectively, T) production data is shown 
in Fig. 6 (respectively. Fig. 7). An excellent agreement is observed on a very wide range 
of center-of-mass energies, spanning from ^/s = 16.8 GeV to ^/s = 7 TeV, on the full 
Xp domain covered experimentally, and for both J/ip and T production. These results 
indicate that the parametrization (3.11) can be safely used in (3.6) in order to compute 
heavy-quarkonium nuclear suppression. 

Independently of the present work, this parametrization can prove useful in future 
phenomenological studies for which a data-driven knowledge of the Xp single differential 
cross sections is necessary. 

^^These data fits have been made much easier thanks to the Quarkonii database which can be found 
at http: //hepdata. cedar. ac .uk/review/quarkonii/. 

We were unable to estimate the experimental errors when retrieving NA3 data from [12] . Therefore we 
do not quote the uncertainties on the exponent n in this case. 
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Figure 6. Comparison between i /ip production data in p-p, 7r^-p and p-A collisions and the fit 
(3.11) (solid red line). The values of n obtained from the fit are indicated in each panel and in 
Table 1. Data are taken from [12, 31-34]. 




-0.2 0.2 0.4 0.6 0.8 -0,3 -0.2 -0.1 0.1 0.2 0.3 0.02 0.04 O.OB 0.08 0.1 



Xp Xp Xp 

Figure 7. Comparison between T production data in p-p and p-d collisions and the fit (3.11) 
(solid red line). The values of n obtained from the fit are indicated in each panel and in Table 2. 
Data are taken from [35-37]. 



3.3 Transport coefficient and nuclear broadening 

The amount of medium-induced gluon radiation, and hence the strength of ij: suppression 
in p-A cohisions, is controlled by the nuclear broadening Ag^ = in Eq. (1.2). For a 
path length L crossed in the target, the broadening reads 

l\ = qL, (3.12) 

where the transport coefficient q in the target nucleus is related to the gluon distribution 
G{x) in a target nucleon as [38]""^^ 

i='^^0^,.G^.AL.)■^'-^,^G^.). (3.13) 

In the latter equation p is the target nucleon number density, and the scaling violations in 
the running of as and in the evolution of the gluon density are neglected since qL <1 GeV^ . 

^*We assume that the relevant transport coefficient in our approach is that of an energetic gluon, lience 
the factor Nc in the numerator of (3.13). 
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The typical value of x at which xG{x) should be evaluated in the r.h.s. of (3.13) is 
discussed in Appendix A. The result depends on whether the hard subprocess is incoherent, 
^hard ^ -^j or Coherent, thard ^ in the nucleus. Assuming 2 — 1 subprocess kinematics 
we estimate thard ~ E/M"^ ~ xiE^/ {xiX2s) ^ l/(2mpX2), so that the incoherent and 
coherent limits correspond respectively to X2 ^ xq and X2 ^ xq, with xq = xq{L) = 
{2nipL)^'^ . In the incoherent case xG{x) should be evaluated at x ~ xq [38], whereas in 
the coherent case x ~ X2 (see Appendix A), i.e., x = min(xo,X2). 

Using the power-law behavior xG{x) ~ x"'^ '^ suggested by small-x (x < 10~^) fits to 
HERA data [39] , we can thus extract the x and p dependence of q by writing 

) ; X = min(xo,X2) ; xq = 7: (3.14) 

X J 2mpL 

where po is in principle an arbitrary constant density, and % = q{x = 10"^, /3 = po). 

In Ref. [19] we used the hard sphere (HS) approximation^^ and thus the same uniform 
density p = p-^g = [(4/3)7rrQ]~^ ~ 0.17 fm~^ for all nuclei. Within this approximation, 
q = q{x) with % = q{x = W~^,p = p^g), the average of L = J dz is found to be 
Ljjg = 3-Ryi/2, and the broadening is directly obtained from (3.12). In the present study we 
will use more realistic (non-uniform) nuclear density profiles p{r) extracted from electron- 
proton scattering experiments [40]. In order to formally recover the situation considered 
in Ref. [19] when p{r) is constant, we thus choose po = Pjjg in (3.14). 

When p{r) is not constant, the parameter L = J dz entering (3.12) is badly defined, 
but the broadening is still well-defined, since it is proportional to J pdz rather than to 
J dz. Using (3.14) we can write 

il = jdil = j qdz = q{x) ■ L,s ] L,s = j pdz. (3.15) 

The effective path length L^s is mathematically well-defined and can be related to the 
number A'^part of nucleons participating to the broadening of the fast color charge. Using 
dA^part = pcrdz, where a is interpreted as the cross section for having non-zero broadening 
in parton-nucleon scattering, we obtain^'^ 

Leff - Lp = ^^PHi^ . (3.16) 
Pocr 

For minimum bias p-A collisions, the average of Apart in the events with J/^ produc- 
tion can be calculated within Glauber theory and reads 

(Apart),/, = 1 + ^^^^^ I d%TA{h)\ (3.17) 

where we used the normalization j d^r p{r) = j d?h TA{b) = A. The effective path length 
becomes 

Les -Lp = ^^2^ / TAibf . (3.18) 

^^This approximation is defined by p(r) = [A/V) e(i?A - r), with Ra = ro A^''^ (ro = 1.12 fm). 

^"The integration constant in (3.16) follows from the fact that for a proton target, A^part ~ 1, see (3.17). 
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Nucleus 


P 


Be 


C 


Ca 


Fe 


Cu 


W 


Pt 


Au 


Pb 


Atomic mass 


1 


9 


12 


40 


56 


63 


184 


196 


197 


208 


Lcff (fm) 


1.5 


3.24 


3.94 


5.69 


6.62 


6.67 


9.35 


10.85 


10.21 


10.11 


^Hs (fm) 


1.68 


3.49 


3.85 


5.75 


6.43 


6.68 


9.56 


9.76 


9.78 


9.95 



Table 3. Values of icft for various nuclei, using (3.19) and realistic nuclear densities [40], or within 
the hard sphere approximation. 

As can be seen from Eq. (3.18), the additional (effective) path length in a nucleus with 
respect to a proton is independent of cr, and can thus be uniquely determined knowing the 
nuclear density profile. For the effective length in the proton, we take Lp = 3i?p/2 = 1.5 fm, 
using Rp = 1 fm for a generic proton length scale. In summary, L^g is given by 

Leff = 1.5 fm + / d% TA{hf . (3.19) 



The values of Leg obtained from (3.19) using realistic nuclear density profiles [40] are listed 
in Table 3 for various nuclei, and compared to the values L-^^ = 2)Ra/2 corresponding to 
the hard sphere approximation previously used in Ref. [19]. 

Using (3.19) in (3.15) fully determines the nuclear broadening and hence the induced 
gluon spectrum (1.2). The transport coefficient go = (li^ = 10~^,/9 = po) is the only free 
parameter of the model. 

3.4 Energy loss probability distribution 

In this section we discuss the quenching weight V{e,E) entering (3.6). A well-known 
procedure (used for instance in the case of large jet-quenching in A-A collisions [7]) to 
construct a normalized V{e,E) from the single gluon emission spectrum dJ/do;, consists 
in assuming independent emissions of soft gluons. In this so-called Poisson approximation. 



oo n 
n=0 U=l 



du 




X 6\e-y]uji] exp{ - I ^ [> . (3.20) 



Let us mention that ^(e, E) is a solution of the equation 

&P{e,E) p dl 



dL 



£ duj[V{e-u:,E)-V{e,E)]^^, (3.21) 



where by convention V{£ < 0,E) = 0, and L stands for any parameter entering the 
expression of dl/duj. If L is the medium length crossed by the fast charge, (3.21) is 
formally identical to the kinetic equation used by Landau to study ionization losses in 
normal matter [41]. An important feature of the Poisson approximation (3.20) is that not 
only each uji, but also the accumulated loss '^uJi, is supposed soft as compared to E. In 



^'^We checked that varying the proton effective length in the range Lp = 1.3-1.7 fm only marginally 
affects our results. 
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other words, the "energy degradation" of the fast particle during the multiple emission 
process is neglected. 

Whether the Poisson approximation is appropriate or not obviously depends on each 
energy loss process and on the specific properties of dl/duj, in particular on the multiplicity 

of radiated gluons with energy ~ O (u). In the present case, dl/duj given in (2.17) is 
a scaling function of uj/Cjj^, where (Da = (^±a/-^^±) -E,^^ and A'^ is estimated as 

~ r du'^r^oj^^ a. In f 1 + 4 V (3-22) 

When uj <^ lDa, becomes potentially large, In ((^a/'^) ~ C)(l)i questioning the as- 
sumption of independent multiple emissions. When lo ^ cja, is small, N^^ < O (a^) <C 
1. However in this region each emitted gluon carries away a fixed fraction ~ i±A/M± of 
the energy E, and the fast particle energy degradation cannot be neglected. Thus, in our 
context the Poisson approximation proves fishy. 

A simple way to deal with this problem is to supplement (3.20) with the condition that 
the energy e is carried away by a single gluon, i.e., 6 {e — X] Wj) — )■ n 5 {e — ujj) in (3.20). 
This yields the (normalized) quenching weight 

P(e.E)=-e.^[-l <i._| = -exp|-y^ d.-| . (3,23) 

The latter is simply interpreted as the product between the "probability" d//de to radiate 
a gluon with uij = £ and the probability (given by the exponential Sudakov factor) to have 
no extra radiation carrying oj^ ^ e. In our context, the expression (3.23) of the quenching 
weight is better founded than the Poisson expression (3.20), and will be used in (3.6). 

3.5 Other nuclear effects 

Besides energy loss effects, other mechanisms might affect ■0 suppression in nuclei. In this 
section, the role of nuclear absorption and saturation on ip suppression in p-A collisions is 
discussed. 

3.5.1 Nuclear absorption 

At small i{j energy in the nucleus rest frame, the hadronization time = ■ [E/M) 
(where is the proper hadronization time) becomes comparable to the typical nuclear 
size, < L. Consequently, ^ states are produced on average within the target nucleus 
and might suffer inelastic interaction with nuclear matter, the so-called nuclear absorption 
process. From (3.3), this should be the case at low proton beam energy E^ [i.e. at low 
^/s ~ y^2mpEp) or at small values of x^. The J/ip suppression in p-A collisions at the 

^^For the present discussion we neglect the second term of the spectrum (2.17). 

^^In the case of large jet-quenching [7], a large multiplicity at small uj is compensated by a factor 
tt/L <IC 1, resulting in a small gluon occupation number ~ (ti/L) N^, <^ 1, thus supporting the assumption 
of independent emissions. In our case the spectrum (2.17) arises from large formation times tf 3> L, and 
the estimates of gluon occupation number and gluon multiplicity coincide. 

^"^This problem was previously addressed in Ref. [42] in the context of large p± jet-quenching. 
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SPS {e.g. by the NA60 experiment at Ep = 158 and 450 GeV [43]) has m particular often 
been attributed to nuclear absorption effects. 

Nuclear absorption effects are not included in this analysis for two reasons. First of 
all, the strength of nuclear absorption strongly depends on the (effective) absorption cross 
section af^^ which is poorly constrained from data [44, 45]. Moreover, when ~ L, the 
hierarchy (2.18) upon which the medium-induced spectrum (2.17) relies is no longer valid 
and hence the use of the latter becomes dubious. In this paper we will therefore focus on 
the region ^ L (the region of validity of (2.17)), where the effect of nuclear absorption 
is irrelevant. 

In the calculations presented in Section 4, we shall indicate by an arrow the value of 
^cnt ycnt^^ defined as E{x'^"^)/M x t^, = L, below which nuclear absorption might play 
a role. For the numerical values of x^""'*, the J/?/' hadronization time is given by the mass 
splitting between IS and 2S states, tj/^ = (M^/ — My^)~^ ~ (0.6GeV)~^ ~ 0.3 fm (a 
similar estimate is obtained in the T channel). Note that x"'* becomes negative at large 
collision energy, in which case the model should not only apply at large positive Xp but 
also down to Xp < 0. 

3.5.2 Saturation and nuclear PDF effects 

At small values of x, partons inside the nucleus wavefunction start to overlap, leading to 
the phenomenon of saturation (see for instance [46] for a review). Although saturation 
effects should also occur in a proton, they are expected to scale roughly like the nucleus 
transverse density, V/S ^ A^^^ , therefore being stronger in large nuclei at a fixed value of 
X. As a consequence, the ip normalized yield in p-A collisions is likely to be suppressed 
with respect to that in p-p collisions - independently of the energy loss effects discussed 
above - either at large Xp and/or at high energies (RHIC, LHC) where small values of x 
are probed in the target nucleus. 

The effects of (gluon) saturation on suppression in p-A and A-A collisions have 
been addressed by many authors, see e.g. [47, 48]. In the present paper, we shall implement 
the physics of saturation following the work of Fujii, Gelis and Venugopalan [48], where 
J/ip suppression has been computed within the Color Glass Condensate assuming 2 — )• 1 
kinematics for the production process. The nuclear suppression is a scaling function of the 
saturation scale Qs and can be simply parametrised as [48] 



with b = 2.65 GeV^ and a = 0.417. Unfortunately, no equivalent parametrization has been 
given in the T channel. We will assume in the present approach that saturation effects 
on heavy-quarkonium production are a scaling function of Qs/M_i.'^^ Therefore, the T 
suppression due to saturation reads 



^^Ideally this ansatz should be checked numerically. Indeed, it can only be approximate, since running 
coupling effects will explicitly spoil this scaling hypothesis. 




(3.24) 



SliQ,) = SriQ.xMi^yMl). 



(3.25) 
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In order to make reliable predictions at RHIC and LHC, the i/ip and T nuclear production 
ratio is determined assuming energy loss effects, B^'}^^^ from Eq. (3.6), with and without 
saturation effects, 

(i) ^pA = ^pa"*"*" ) 

(ii) Rlt = Rft'' X 5a/5p . 

The saturation scale appearing in (3.24) is closely related to the transport coefficient 
q given by (3.13), namely [28, 49] 

Ql{x,L) = qL. (3.26) 

In other words is nothing but the transverse momentum broadening discussed in Sec- 
tion 3.3, see Eqs. (3.12) and (3.15). The inclusion of saturation effects thus does not 
require any additional parameter once the parametrization (3.14) for q is employed and go 
is determined. 

Let us mention that go is related to the saturation momentum in a proton at x = 10"^, 
QsO- We have, from (3.14) and (3.26), 

Qso = Ql{x2 = 10-2, Lp) = go ^p - 0.1 GeV2 , (3.27) 

with go in GeV^/fm. Comparing the value of go obtained in our model from a fit to the 
E866 J/'f/' nuclear suppression data and the current estimates of QIq obtained from a fit 
to small-X2 DIS data [50] should provide a non-trivial (though not conclusive) test of our 
model. 

Another, earlier approach in order to model the modifications of parton densities in 
nuclei is the use of leading-twist nuclear PDF (nPDF) which have been determined from 
global fit analyses of e-A DIS or p-A Drell-Yan data for more than a decade (see e.g. [51]). 
In this framework, ip production in p-A collisions is proportional to the gluon distribution 
in the nucleus G^{x2,M^). Therefore, '0 suppression can be modelled as 

(iii) i?^™^ = X Gf °^(X2, M^)/Gp(x2, M^) . 

The predictions to be discussed in the next section will be performed assuming energy 
loss effects, supplemented with predictions including saturation effects at RHIC and LHC 
energies where these are expected to play a role. For completeness, we will also critically 
compare in Section 4.7 these results with those obtained using nPDF. 

4 Phenomenology 

After the description of the energy loss model, the phenomenology of "0 suppression in 
hadron-nucleus collisions is investigated in this section. In the practical applications, 
we take Aq^^j^ = 0.25 GeV, = 1 GeV in the transverse mass Af_L = -|- p\, and 

M = 3 GeV (M = 9 GeV) for the mass of a compact cc {hh) pair. As we can easily 
verify a posteriori, the typical scale entering the running coupling constant is not too large, 
qL ~ IGeV^, which justifies the assumption of a frozen coupling, = 1/2, at such 
semi-hard scales. 
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4.1 Fitting procedure 

The only parameter of the model, the transport coefficient qo, is determined by fitting the 
J/'0 suppression measured by E866 [9] in p-W over p-Be collisions {^/s = 38.7 GeV). This 
choice is motivated by the fact that the E866 measurements are the most precise performed 
so far and cover a wide range in x^. We choose to perform the fit in the [0.2-0.8] x^- 
range for the following reasons: at x^ < 0.2 J/ip suppression might be affected by nuclear 
absorption (see Section 3.5.1) while at > 0.8, we expect quark-induced subprocesses 
to come into play, possibly modifying the overall normalization of the medium-induced 
spectrum (1.2). Note also that this x^^-range at E866 energy corresponds to values of 
> 10~^ for which saturation effects are expected to be small, of the order of 5% at most 
on the W/Be ratio. 

The fit gives = 0.075 it 0.005 GeV^/fm, where the quoted uncertainty is determined 
from the minimization procedure. A systematic uncertainty on the value of % can be 
roughly estimated by restricting the x^^-range used for the fit to the interval [0.3-0.7]; we 
found that it would increase the value of to % ~ 0.087 GeVVfm. 

The result of the fit is shown in Fig. 8 where excellent agreement is observed in the 
whole fit range. Note however that for x^ ^0.1, nuclear absorption is expected to play a 
role; see the vertical arrow at x^'''* ~ 0.07 below which the J/ip formation time becomes 
smaller than the size of the target tungsten nucleus (see Section 3.5.1). 
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Figure 8. E866 J/tp suppression data [9] in p-W collisions compared to the energy loss model. 

It is worth mentioning that the fitted transport coefficient, % = 0.07-0.09 GeV^/fm, 
would correspond to the saturation scale in a proton Ql{x = 10~^) = 0.11-0.14 GeV^ using 
(3.27), which is consistent with (yet slightly smaller than) estimates based from fits to F2 
DIS data [50] . Note that the saturation scale in large nuclei and at smaller x considerably 
exceeds that in a proton, yielding ~1 GeV^, where the use of perturbative techniques 
is commonly assumed to be legitimate. 
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4.2 Scaling properties of heavy-quarkonium suppression 

Before comparing the model predictions to the other available data, we discuss in this 
section the expected scaling properties of ij: suppression in the present model. 

Let us first mention that the nuclear dependence of quarkonium suppression is often 
parametrised as a power law, 

= ^ A-\ (4.1) 

where a is assumed to be independent of A. The power law is empirical. It can be inferred 
in the Glauber picture of '0 absorption in the nucleus, S'^^^ ~ exp (—est • A^/'^), and using 
the approximation A^/'^ ~ log A, which is accurate to the 10% level for 5 < ^ < 200. 
However, the Glauber picture of nuclear absorption is expected to hold when the ip energy 
E in the nucleus rest frame is small enough, see Section 3.5.1. The heuristic law (4.1) has 
no reason to be valid at high E where the compact color octet QQ pair crosses the nucleus 
and hadronizes far beyond. 

We checked that the i/ip suppression expected in our model does not follow the pa- 
rametrization (4.1). To illustrate this, the typical values of a are found to vary by up to 
10% depending on whether i/if) suppression in p~W collisions (in the E866 kinematics) is 
normalized either to p-p or p-Be collisions. Clearly the attenuation factor R^a should be 
preferred to the effective power a when discussing nuclear suppression. We thus focus on 
i?pA rather than on a throughout our study, and now discuss its scaling properties. 

In the energy loss model of Gavin and Milana [20], quarkonium suppression exhibits 
an approximate Xp scaling, i.e., i?pA is a function of Xp but independent of ^/s. Indeed, 
assuming that the shape of dcjpp/dx^ is independent of -y/s, and considering the limit 
Xp » M^l^ where E XpEp, we obtain from (3.7) 



In the present approach, however, the approximate Xp scaling of quarkonium suppression 
is broken for several reasons: 

1. The transport coefficient q, and therefore the function Jioss in (4-2), depends expli- 
citly on X2 at small X2 < xq, see (3.14). As we shall see, this effect is particularly 
important at LHC energies; 

2. As discussed in Section 3.2, the slope of the p-p production cross section does depend 
on y/s (see Tables 1 and 2); 

3. Finally, the saturation (or nPDF) effects also scale with X2 yet this effect is actually 
rather small. 

In order to illustrate this, J/V' suppression has been computed in p-W collisions as a 



function of Xp at NA3 = 19.4 GeV), E866 {^/s = 38.7 GeV), RHIC {^/^ = 200 GeV) 




(4.2) 
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for various assumptions regarding the transport coefficient q and the value of the exponent n. See 
text for details. 

and LHC (y^ = 5 TeV) energies in Fig. 9, without saturation effects but under various as- 
sumptions. In the left panel, the transport coefficient is frozen, q{x) = % = 0.075 GeV^/fm, 
and the exponent of the p-p cross section is fixed to n = 4.5. With no surprise the scaling 
is observed to a very good accuracy, except at the lowest ^/s for which the approximation 
E XpEp is no longer valid. When taking explicitly into account the x dependence of q{x) 
but keeping a fixed exponent n = 4.5 (central panel), the Xp scaling is strongly violated 
at LHC energy, but still approximately verified from NA3 up to RHIC energies. Finally, 
the deviations from Xp scaling are even more pronounced (right panel) when considering 
the actual exponents n extracted at each y/s in Section 3.2, with a stronger suppression at 
RHIC (n = 8.3) and LHC (n = 32.3). 

Note that at LHC, the variation of R^Jp^ with Xp is extremely fast at very small Xp. 
This strong dependence comes from the small-x behavior of the transport coefficient q{x), 
Eq. (3.14), together with the fast variation oi x = x^ with Xp at \xp\ < 10~^. As we shall 
see in Section 4.5 the variation of R^Jp^ with y ^hiXp is naturally much smoother. 

In order to check experimentally whether J/?/' suppression scales with Xp, it would be 
crucial to measure J/V' production in p-A collisions at RHIC and LHC at large x^, say 
Xp > 0.1, which is out of reach with the present apparatus. Such measurements could in 
particular shed light on the x dependence of the transport coefficient q{x). 

4.3 Predictions and comparison to J/V' data 

Once go is determined from the fitting procedure described in Section 4.1, the Xp depen- 
dence of the i/ip quenching factor i?p^ can be predicted in any target nucleus and at any 
center-of-mass energy for which the absolute p-p differential cross section dcjpp/dx^ has 
been measured. In this section we systematically compare the model predictions with all 
available data. 

4.3.1 E866, NA3, E537, NA60, HERA-B 

Let us start with the comparison oi i/tjj suppression expected in an iron target and the 
E866 data for iipc/Bei ^-C-, taken at the same energy as the fitted ratio -Rw/Bo- The excellent 
agreement reported in Fig. 10 fully supports the atomic mass dependence of the model. 
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This is at variance with the calculations by Gavin and Milana [20] which overestimated 
the ratio -Rpe/Bej that time measured by E772.^^ It is therefore a hint that the L- 
dependence expected here, AE oc ^/L (see (1-3)), is probably more appropriate than the 
ad hoc assumption of Ref. [20], AE oc L. 

Data taken at lower or smaller Xp are also compared to the model. As can be seen 
in Fig. 11 the agreement with NA3 p-A and 7r~-A data is excellent, both in shape and 
magnitude, over a very wide range in Xp. It is also remarkable that the model is able to 
reproduce the different magnitude of suppression in p-A and 7r~-A collisions reported by 
NA3 [12]. This difference cannot be understood within nuclear absorption models, where 
nuclear suppression is a purely final state effect, thus independent of the projectile type. 
It cannot either be explained by nPDF effects, unless the nPDF to proton PDF ratios 
for valence quarks and for gluons, probed respectively in tt~-A and p-A collisions, prove 
completely different. In our picture, the smaller J/ip suppression in 7r~-A collisions 
naturally arises from the flatter differential cross section, jiTrp = 1.4 vs. npp = 4.3 at 
-y/s = 19.4 GeV, see Table 1. Although no prediction of the exponent n is made in our 
model, it is clear that this feature can be explained from the larger slope, at large x, of 
the gluon PDF in a proton, xG{x) ~ (1 — x)^ [54], when compared to that of a valence 
antiquark PDF in a pion, xq{x) ~ [1 — x) [55]. In this respect, let us mention that our 
assumption of an incoming gluon in quarkonium hadroproduction (see the Introduction) 
does not hold for NA3 pion-nucleus collisions, where subprocesses with an incoming valence 
antiquark dominate. In spite of this, a very good agreement between the model and the 
NA3 7r~-A data is found, suggesting a mild dependence of the energy loss on the incoming 
parton type. A similar remark applies to the case of the 7r~~A E537 data discussed below. 

The E537 experiment also reported on measurements oi J/tp production in vr" induced 
collisions on various nuclear targets (Be, Cu, W) at ^/s = 15.3 GeV [56]. Our results^^ are 
found in reasonable agreement with the measured ratio -Rw/Be (Fig- 12, left); the slight 
underestimation of the suppression by the model might be attributed to nuclear absorption. 
Indeed, at this energy x^^^^ ~ 0.7, and all E537 data lie in the Xp < x^^'* domain. This 
might also explain the (more pronounced) difference between the observed and predicted 
magnitudes of the ratio -Rw/Cu (Fig. 12, right). 

In Figs. 13 and 14 we compare our predictions with NA60 [43]^^ and HERA-B [32] p-A 
measurements. Although the center-of-mass energy is larger than those of NA3 and E537, 
the typical .l/ip energy range covered by NA60 and HERA-B is actually lower because of 
the smaller Xp values probed by these experiments. As a consequence, 3/iJj suppression can 
be affected more strongly by nuclear absorption effects, as can be inferred by the position 
of the x"'* arrows in Figs. 13 and 14, below which hadronization typically takes places 
inside the nuclear target. ^'^ 

^®We do not show the agreement between our model predictions and J/ip E772 data [52] since those 
measurements were superseded by E866 [9]. 

■^^On top of this, nPDF effects in the NA3 kinematical domain, X2 ~ 0.1-0.2, are known to be small for 
both valence quarks and gluons, see for instance the discussion in [53]. 

^^Lacking Tv~-p data at E537 energy, we choose the exponent n = 1.4 (see Table 1). 

^^Lacking p-p data at NA60 energies, we choose the exponent n — 4.3 (see Table 1). 

■^"in the left panel of Fig. 13, the arrow is not visible as a;^'''* > 0.4. 
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Figure 10. E866 J/-0 suppression data [9] in p-Fe collisions compared to the energy loss model. 
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Figure 11. NA3 J/i/' suppression data [12] in p-A and 7r~-A collisions compared to the energy 
loss model. 

Nevertheless, the model predictions prove in very good agreement with data. In par- 
ticular, the enhancement of J/V' production observed by HERA-B at very negative Xp, 
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Figure 12. E537 J/'0 suppression data [56] in tt -A collisions compared to the energy loss model. 
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Figure 13. NA60 J/i/j suppression data [43] in p-A collisions compared to the energy loss model. 



Xp < —0.2 (see Fig. 14) is well reproduced by the model. (The origin of RpA > 1 can be 
simply understood from the positive slope of dcjpp/dx^ in the target fragmentation region, 
Xp < 0, see the HERA-B data in Fig. 6.) There is however room for J/i/^ absorption with a 
cross section of a few millibarns, as suggested by the slight overprediction of RpA at NA60 
precisely in the region where J/ip absorption can no longer be neglected. 

4.3.2 PHENIX 

The predictions in d-Au collisions at RHIC, y/s = 200 GeV, are shown in Fig. 15 in 
comparison with PHENIX data [13], with (dashed line) and without (solid line) saturation 
effects. The energy loss model is able to reproduce nicely the J/^ suppression at all 
rapidities. Note that in several phenomenological analyses, the suppression observed in 
the most forward rapidity bins has often been attributed to gluon saturation effects or to 
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Figure 14. HERA-B suppression data [32] in p-A collisions compared to the energy loss 
model. 

strong small-x shadowing in the nuclear PDF (see e.g. [47, 57]). Here, the sole energy loss 
effects might be responsible for the observed suppression, although saturation might play 
a role as well. As a matter of fact, the agreement is better when saturation is included. 
Remarkably, an excellent agreement is also observed in some negative y bins, for which 
nuclear absorption might also play a role (at least for y < y'^^'* = —1.1). We shall discuss 
further these data in Section 4.7 when comparing to the predictions including nPDF effects. 

The J/V' suppression has also been measured by PHENIX for various d~Au centrality 
classes [13] and more recently as a function of its transverse momentum [58]. Discussing 
these data would go beyond the scope of the present article and is left for future work. 

4.4 Predictions and comparison to T data 

The above comparison between J/ip suppression data and our model predictions supports 
both the medium length and energy dependence of the model. The mass dependence 
of heavy-quarkonium suppression can be studied by investigating the suppression of T 
states in p-A collisions. Unfortunately the data are rather scarce; to our knowledge, 
the measurements have only been performed by E772 at Fermilab [59] and PHENIX and 
STAR [36, 60] at RHIC. 

The E772 data are shown in Fig. 16 for various nuclear targets (Ca, Fe, W) and in 
comparison to the model predictions. A rather good agreement between data and theory 
is found for Xp > x™'*, although smaller experimental uncertainties would be necessary to 
further check the M dependence of the model. At small Xp < x™*, the measured nuclear 
production ratio lies much below our predictions, probably too much to be accommo- 
dated by T nuclear absorption. However, let us mention that the E772 measurements at 
low Xp might be affected by uncorrected acceptance effects due to the correlation between 
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Figure 15. PHENIX J/V-' suppression data [13] in d-Au collisions compared to the energy loss 
model, with (dashed line) and without (solid line) saturation effects. 

Xp and (see the discussion in [9]); it is therefore difficult to draw any conclusion from 
the significant disagreement observed at negative Xp. 
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Figure 16. E772 T suppression data [59] in p-A collisions compared to the energy loss model. 

These T data nevertheless allow the mass dependence of the energy loss to be con- 
strained. In their paper [20], Gavin and Milana assumed that the mean energy loss scales 
as /\E oc M~", and considered explicitely the cases n = 2 ("power suppressed") and n = 0. 
From the comparison of their calculations with E772 data, these authors concluded that 
neither of these two choices were satisfactory: assuming /\E oc led to too little T 

attenuation while a too strong suppression was predicted with the hypothesis AE oc M^. 
It is therefore interesting to note that the scaling AE oc predicted in [15] and used 
here (see (1-3)) supports this empirical observation, as the agreement in Fig. 16 indicates. 

The predictions at RHIC are shown in Fig. 17. As expected, the suppression is less 
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Figure 17. PHENIX (|y| = 1.7) and STAR {y = 0) T suppression data [36, 60] in d-Au collisions 
compared to the energy loss model. 

pronounced than for .l/ip production, compare to Fig. 15. Since saturation effects are very 
small in the T channel, the predictions including saturation or not are virtually indistin- 
guishable. The PHENIX and STAR experiments reported on preliminary measurements 
of T suppression^^ in d-Au collisions [36, 60], see Fig. 17. Hopefully more precise data will 
soon allow for clarifying the strength of T suppression in cold nuclear matter. 

4.5 LHC predictions 

We discuss in this section the J/ip and T suppression expected in p-Pb collisions at the 
LHC {y/s = 5 TeV). In Fig. 18 we show the RpPh ratios for both states as a function of 
the rapidity in the center-of-mass frame. Interestingly, J/ijj production is significantly 
suppressed at large positive rapidity, e.g. -Rpp^ — 0.7-0.8 at y = 1 and down to i?ppb ^ 0.5 
y ^ 4- Because of the high center-of-mass energy of the collision at the LHC, saturation 
effects in the J/ip channel are significant: in addition to energy loss, the suppression due to 
saturation ranges from sj^'^ ~ 0.9 in the most negative rapidity bins down to S'j^'^ ~ 0.65 
at y = 5. In the target fragmentation region (y < 0), the suppression is moderate (~ 10- 
20%) while a possible J/ip enhancement might be seen at very backward rapidities, y < —5. 
Predictions using EPS09 [61] and DSSZ [62] nPDF sets are also discussed in Section 4.7. 

In the T channel, the suppression is more moderate because of the mass dependence 
of energy loss,^^ AE oc M_l"\ e.g. R^p^ ~ 0.85 at y = 3. At the LHC the saturation 
effects in the T channel prove quite small, although more pronounced than at RHIC. As 
can be seen from the arrow in Fig. 18, J/ip and T hadronization should take place outside 

^^The PHENIX and STAR data correspond to the suppression of T(1S)+T(2S)+T(3S) states. 
^^Note that in p-Pb colhsions at the LHC, the laboratory frame is shifted by Ay ~ 0.47 with respect to 
the center-of-mass frame. 

Another effect, yet rather marginal, comes from the flatter distributions in p-p colhsions (see Table 2 
in Section 3.2). 
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Figure 18. J/i/' and T suppression expected in p~Pb collisions at the LHC, with and without 
saturation effects (see legend). 

the nuclear medium for y > ~ —5; nuclear absorption should thus play little or no 
role at the LHC. 

These predictions can be compared to the future measurements during the p-Pb run 
scheduled at the LHC in January 2013. In order to test the model, the nuclear dependence 
of ifj production should ideally be measured for various rapidity bins and on a rather broad 
range, which hopefully should be possible with the ALICE or LHCb experiments.^^ 

4.6 E906 predictions 

The E906 "SeaQuest" collaboration [63] aims at measuring Drell-Yan production in p-p 
and p-A collisions at Ep = 120 GeV {^/s = 15 GeV) at Fermilab. Although the first goal of 
this experiment is to study the sea quark asymmetry in the nucleon, it will also be able to 
measure the nuclear dependence of both Drell-Yan and J/ip production on various nuclear 
targets and on a wide range in . In this section we present our model predictions on J/ip 
suppression in p-A collisions at E906 energy, to be compared to the measurements that 
might already be available in 2013. 

In Fig. 19 we plot the predictions in p-Fe (left) and p-W (right) collisions. The 
suppression is very pronounced especially at large Xp , for which however the J/ip production 
cross section should be extremely small. 

^*We recall that the calculations are made using — 1 GeV in the calculation of the transverse mass. 
Therefore our predictions on J/tp suppression should be adapted for the CMS acceptance which requires a 
transverse momentum cut, Pj^ > 6 GeV, in the J/ip channel. 

^^Lacking p-p data at this energy, we choose the exponent n = 4 to be consistent with the systematics 
discussed in Section 3.2. 
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Figure 19. i/ij] suppression in p-Fc and p-W collisions in the E906 kinematics. 
4.7 Comparing predictions using saturation vs. nPDF 

For completeness, we compare in this section the former results on J/'i/; suppression at 
RHIC and LHC obtained in the "energy loss + saturation" model with the predictions 
using the EPS09 [61] and DSSZ [62] nPDF leading-order sets instead of saturation. Unlike 
saturation effects, nPDF corrections should be valid (and possibly non-negligible) even at 
not too small values of x^, and in particular at E866 energy. Therefore, the transport 
coefficient qo using each of the two nPDF sets has been consistently refitted to E866 data. 
The corresponding values, used for the RHIC and LHC predictions with these two sets, are 
indicated in Table 4. 

The comparison is shown in Fig. 20 at RHIC. The predictions using saturation, EPS09 
and DSSZ somehow differ in the rapidity dependence of -R^Au- "^^^ DSSZ modifications are 
rather small, leading to a suppression similar to the one assuming energy loss effects only. 
On the contrary, the EPS09 set exhibits larger modifications to the gluon nPDF (and in 
particular a slightly faster variation in this x domain) increasing the slope of -Rj/^ versus y. 
At mid- and forward rapidity, the various predictions are similar; in particular all of them 
reproduce nicely the data at y = 1~2 with a slightly better description with saturation or 
using the EPS09 set as compared to DSSZ (yet this is not statistically significant). On 
the other hand, the predictions are different at backward rapidity. The best agreement is 
obtained assuming saturation effects (in addition to parton energy loss) or DSSZ nPDF 
instead of EPS09. This observation obviously depends on the present energy loss model, 
preventing us from drawing a firmer conclusion. 

In order to analyze a bit more quantitatively these results, the values of x^/ndf ob- 
tained for the E866 and PHENIX data sets are quoted in Table 4. As can be seen, the E866 
data do not allow one to differentiate the various energy loss predictions including (or not) 
EPS09/DSSZ nPDF corrections. On the contrary, the agreement at RHIC is considerably 
better when energy loss is supplemented by saturation effects (x^/ndf = 0.3) rather than 
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Figure 20. J/tp suppression predicted in d-Au collisions at RHIC in the energy loss model, for 
various assumptions regarding the nuclear modifications of gluon distributions in nuclei. PHENIX 
data are from [9]. 
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Figure 21. Left: J/ip suppression predicted in p-Pb collisions at the LHC in the energy loss model, 
for various assumptions regarding the nuclear modifications of gluon distributions in nuclei. Right: 
same when normalized to its expected suppression at mid-rapidity, RpPh{y) / RpPhiy ~ 0). 



by nPDF (x^pg^g/ndf = 2.7, x^ssz/'^'^^ ~ ^■^): as mentioned above. For completeness we 
also quote the values of x'^/ndf assuming no energy loss but only saturation"^^ or nPDF 
corrections. As can be seen from Table 4, saturation without energy loss still gives a fair 
description of PHENIX data, x^/ndf = 2.7 (as well as EPS09 and DSSZ to a lesser extent), 
although saturation and nPDF effects alone would totally fail to reproduce E866 data. 

■^^The value of go = 0.075 GeV^/fm quoted in Table 4 is here to determine the saturation scale. 
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Table 4. and x^/ndf of E866 and PHENIX data, with (upper rows) and without (lower) energy 
loss effects, for various assumptions regarding the nuclear modifications of gluon distributions in 
nuclei. 

The predictions at the LHC are shown in Fig. 21 (left). The expected suppression 
with saturation effects or using the EPS09 set prove remarkably similar. On the contrary, 
the nPDF corrections given by DSSZ are tiny in the forward rapidity bins, despite the 
small values of x probed in the Pb nucleus. Although the absolute magnitude of the J/i/^ 
suppression differs depending on the assumption regarding saturation/nPDF effects, the 
rapidity dependence (especially at y > 0) is mostly governed by energy loss effects. This 
could be tested experimentally. In the present model, energy loss effects are moderate at 
mid-rapidity which corresponds to the maximum of the p-p production cross section. As 
a consequence, the expected suppression at y = 0, i?ppb(y = 0), is more sensitive to satu- 
ration/nPDF effects. Moreover, since the rapidity dependence is essentially due to energy 
loss effects, the double ratio -RpPb(y)/-RpPb(y = 0) is rather independent of the strength of 
saturation/nPDF effects. This is illustrated in Fig. 21 (right) where RpPhiu) / RpPhiu = 0) 
is plotted. As can be seen this double ratio proves remarkably similar whether or not 
energy loss is supplemented with nPDF or saturation effects. 

5 Discussion 

The agreement between our model and the p-A data for quarkonium nuclear suppression 
is quite remarkable. With a single free parameter %, both the slope and normalization of 
RpA (or Rp^/RpB and also R-ka) are accurately described, for various collision energies, 
various target nuclei and different masses {<l/tp, T), and over a broad range in Xp (or 
rapidity) . We also stressed that the effect of saturation alone fails in describing J nuclear 
suppression at different collision energies. This strongly supports parton energy loss as a 
dominant effect in p-A quarkonium nuclear suppression, the main conclusion of our study. 
The successful description of the data is mostly due to the (medium-induced) energy loss 
scaling as AE oc E, where E is the energy of the QQ pair in the nucleus rest frame. 
This behaviour arises when the partonic subprocess looks like small angle scattering of 
an asymptotic charge, and thus holds within our assumption of a long-lived, color-octet 
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QQ pair. Our results support the parametric {E, M and L) dependence of the induced 
radiation spectrum (1.2), which is derived from first principles in Section 2. 

These results also give some hint on the mechanism for low p±_ {p± < M) heavy- 
quarkonium hadroproduction. We argued in the Introduction that at large Xp, the octet 
QQ pair should be long-lived in any quarkonium production model, including the Color 
Singlet Model (CSM). The agreement of the energy loss model with the large sup- 
pression data thus cannot distinguish between production models. But we found that the 
agreement extends to small values of Xp (see in particular the comparison to RHIC data 
at y ~ 1-2, corresponding to Xp ~ 0.04-0.1, in Fig. 15), where assuming a long-lived color- 
octet QQ pair becomes inaccurate in the CSM. The CSM mechanism thus seems somewhat 
disfavoured by our results, at least as a dominant contribution to inclusive {i.e., low p± 
and low Xp) J/^ production. The future measurements in p-Pb collisions at the LHC will 
probe small values of Xp {\Xp\ < 0.1), yet in a rather large rapidity interval {\y\ < 5), and 
might thus further clarify the underlying dynamics of heavy-quarkonium production. In 
fact our energy loss explanation of J/V' suppression is consistent with any J ftp production 
model where ihard ^ ^octct, leaving room for gluon radiation with thard ^ tf ^ toctet, see 
(2.18). It was argued in Ref. [30] that a qualitative analysis of the quarkonium production 
data suggests a mechanism for quarkonium hadroproduction, named "Comover Enhance- 
ment Scenario", where color neutralization is realized at the time toctet by a semi-hard 
scattering between the QQ pair and the comoving radiation field of the incoming par- 
ton. It is intriguing that the condition on toctet inferred from global production features, 
^hard ^ ^octct [30], is Consistent with the condition (2.18) necessary to explain nuclear 

suppression from radiative parton energy loss. 

The induced radiation spectrum was derived assuming a given hierarchy between the 
nuclear size L, gluon formation time tf and quarkonium hadronization time t^. Thus, 
as we emphasized several times, the model should in principle be valid only when the 
quarkonium state hadronizes outside the nucleus, i.e., when E is large enough or Xp > x'^p^^. 
It is quite striking that the extrapolation of the model to the region Xp < x^^^^ is either 
consistent with the data (within error bars, see e.g. the NA3 data in Fig 11, HERA-B data 
in Fig. 14 and PHENIX data in Fig. 15), or systematically underestimates quarkonium 
nuclear suppression (NA60 data in Fig. 13). This suggests parton energy loss to play a 
role in a broader domain than expected, the possible additional suppression required at 
Xp < x^"^ being due to nuclear absorption of the fully formed quarkonium state. 

In our study we also assumed quarkonium production in proton-nucleus collisions to 
arise from the splitting of an incoming gluon into an octet QQ pair. This assumption 
becomes inaccurate at very large Xp, where quark-induced processes (such as qq — QQ) 
come into play. Although we expect the parametric dependence of the associated radiation 
spectrum to be unchanged, the overall color factor might be changed in this case. A possibly 
smaller effective color factor at very large Xp might explain the milder J/iJj suppression 
observed by E866 at Xp ^ 0.8 (see Fig. 8) than predicted in our model. However, as already 
mentioned in Section 4.3.1, the very good agreement between the model and the NA3 
7?«on-nucleus data, for which the qq annihilation channel is dominant at all Xp , suggests a 
relatively weak dependence of the energy loss on the incoming parton type. 
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We might envisage refinements of tlie parton energy loss model presented here, such 
as including quarkonium absorption at Xp < x'^^^^ and quark-induced processes, in order 
to extend the domain of validity of our approach. However, we find it more important to 
first confirm the dominant role of parton energy loss in p-A collisions, where gluon-induced 
processes dominate and our model assumptions apply. 

First, the energy loss model can be tested in forthcoming p-Pb collisions at the LHC, 
for which our predictions for the y-dependence of J/ip and T suppression are shown in 
Fig. 18, and in p-A collisions in the E906 fixed-target experiment at Fermilab (Fig. 19). The 
model should as well be confronted to the existing RHIC d-Au data on the p± and centrality 
dependence of J/^ suppression, measured at various rapidities [13, 58]. This requires 
generalizing (3.6) to double differential (in Xp and p±) cross sections, and will be the subject 
of a future work. It will be interesting to check whether the L-dependence of the energy loss 
predicted in (1.3) is consistent with the centrality dependence of the RHIC d-Au data. Our 
study might also help interpreting quantitatively quarkonium measurements performed in 
heavy-ion collisions at RHIC [64, 65] and LHC [66-68] . Indeed, parton energy loss through 
the incoming cold nuclei is expected to combine with hot effects (such as Debye screening 
or final state energy loss in a QGP). The evaluation of J/ip suppression in A-A collisions 
expected from cold parton energy loss alone will be presented in a future study. Finally, 
as discussed in the Introduction other processes than quarkonium production should be 
sensitive to a parametrically similar parton energy loss, such as open charm and light 
hadron production in p-A collisions. This work is also in progress. 
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A X dependence of q 

In Ref. [38], the transport coefficient q was related to the gluon distribution G{x) in a 
target nucleon, see the expression (3.13). The value of x to be used in xG{x) in (3.13) can 
be estimated by considering the kinematics of the rescattering process. 

Following Ref. [38], let us consider the specific case of DIS, where an energetic light 
quark of momentum p is produced and then rescatters with transfer g on a target nucleon 
of momentum P. Working in the target nucleus rest frame, we have P = (mN,C)), with 
TTT-N the nucleon mass. Choosing light-cone coordinates p^ = p^ ±p^ and p = {p~^ ,p~ ,0±) 
along the negative z-direction, the condition for the final quark to be on-shell reads 

[p + qf = {p+ + q+){p- +q-) - ql = ^ p+ + q+ ^ (A.l) 

p 

where we neglected q~ as compared to p" = 2E. 
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parton produced inside the target 

When the hard production time thard ~ E/Q"^ <^ L, or equivalently when the Bjorken 
variable xb = Q"^ / {2m^E) ^ xq = 1/ {2m^L), the parton p is effectively produced inco- 
herently inside the nucleus. This is the situation considered in Ref. [38], which we now 
briefly review. If the rescattering occurs at a distance z from the production point, from 
the uncertainty principle we have ~ 1/z just before the scattering. From the con- 
straint z < L, we obtain ^ 1/L. For p~ large enough Eq. (A.l) gives ~ \p~^\- The 
momentum fraction of the rescattering gluon thus satisfies [38] 

parton produced far before the target 

When thard ^ L <^ Xb ^ xq, the virtual photon splits into a light quark-antiquark pair far 
before the nucleus, and the DIS process is coherent over the whole nucleus. In this case, 
the quark virtuality = ~ Q^, and \p'^\ ~ Q'^/p~~ is not bounded by 1/L any 

longer. From Eq. (A.l) we obtain (using q'j_ <^ Q^) 

X = ~ ~ = XB [XB < Xq) . (A. 3) 

mN mN ?TiN p 

For a generic hard process (for instance in a p-A collision) of coherence length thard ~ 
E/M"^ ~ l/(2mN2;2), the above DIS example supports the following estimate for the value 
of X to be used in Eq. (3.13), 

x = xqQ{x2> XQ) + X2Q{x2<xo)=mxa{xQ,X2) ; Xq = —^—-, (A.4) 

thus specifying the 3;2-dependence of the transport coefficient q. 
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